Open-system dynamics of graph-state entanglement 
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We consider graph states of arbitrary number of particles undergoing generic decoherence. We 
present methods to obtain lower and upper bounds for the system's entanglement in terms of that 
of considerably smaller subsystems. For an important class of noisy channels, namely the Pauli 
maps, these bounds coincide and thus provide the exact analytical expression for the entanglement 
evolution. All the results apply also to (mixed) graph-diagonal states, and hold true for any convex 
entanglement monotone. Since any state can be locally depolarized to some graph-diagonal state, 
our method provides a lower bound for the entanglement decay of any arbitrary state. Finally, this 
formalism also allows for the direct identification of the robustness under size scaling of graph states 
in the presence of decoherence, merely by inspection of their connectivities. 



PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Yz 

Introduction.- Graph states pQ constitute an impor- 
tant class of entangled states with broad-reaching appli- 
cations in quantum information, including measurement- 
based quantum computation [21 [3J, quantum error correc- 
tion U, and secure quantum communication [5]. More- 
over, instances of this family, such as the Greenberger- 
Horne-Zeilinger states, play a crucial role in fundamental 
tests of quantum non-locality . Consequently, a great 
effort has been made both to theoretically understand 
their properties Q] [7] and to create and coherently ma- 
nipulate them experimentally [8^ 

Needless to say, it is crucial to understand the dynam- 
ics of their entanglement in realistic scenarios, where the 
system unavoidably decoheres due to experimental errors 
or to the interaction with its environment. Previous stud- 
ies on the robustness of graph-state entanglement in the 
presence of decoherence observed a disentanglement time 
(or lower bounds thereof) insensitive to the system size 
fTUj . However, the disentanglement time on its own 
is not in general able to provide any faithful assessment 
about the entanglement's robustness, since it can grow 
with the number TV of particles and yet the entangle- 
ment can get closer to zero the faster, the larger N [IT] . 
The full dynamical evolution of entanglement must then 
be studied to draw conclusions on its fragility. Taking the 
latter into account, the entanglement of the linear-cluster 
states, an example of graph states, was shown to be ro- 
bust with the size of the system against the particular 
case of collective dephasing decoherence [T5] . 

The present work provides a general framework for 
the study of the entanglement evolution of graph states 
under decoherence. Our techniques apply to (i) any 
graph, and graph-diagonal, states; (ii) arbitrary kinds 
of noise, individual or collective; and (iii) any convex (bi- 
or multi-partite) entanglement quantifier that does not 
increase under local operations and classical communi- 



cation (LOCC). In the developed formalism we consider 
local measurement protocols to efficiently obtain lower 
and upper bounds for the entanglement of the whole sys- 
tem contained in any given partition in terms of that of a 
considerably smaller subsystem consisting only of those 
qubits lying on the boundary of the partition. No opti- 
mization on the full system's parameter space is required 
throughout. For an important class of noisy channels - 
namely arbitrary Pauli maps, to be defined below - the 
lower and upper bounds coincide, providing thus the ex- 
act entanglement evolution. With the same methods we 
also establish a second family of lower bounds that, de- 
spite less tight, depend only on the connectivity of the 
graph and not on its size. This allows us to assess the ro- 
bustness based on the full dynamics of the entanglement 
and not just its disentanglement time. Our approach can 
also be used to establish lower bounds to the entangle- 
ment behavior of any initial quantum state. 

Graph states.- Consider a mathematical graph 
G(y,f) = {V,£}, composed of a set V, of N vertices 
i € V, and a set £, of edges {i,j} £ £ connecting each 
vertex i to some other j. The associated physical state 
is operationally defined as follows: to each vertex i as- 
sociate a qubit, initialize all N qubits in the product 
state | 5(v)o ) = l£V being 1+,) = (|0<) + |1;))/V2, 
and to all pairs {i,j} of qubits joined by an edge ap- 
ply a maximally-entangling control- Z (CZ) gate, CZij = 
1 0,0^ ) (Oi 0^ I + 1 0, 1 ^ > (0, 1 ^ I + 1 1,0^ ) < 1 ,0^ I - 1 1, ) < 1 j I ■ The 
resulting iV-qubit graph state is 

l G (V,£) )= <S> CZ iMv) )- (1) 

An example of such graph is shown in Fig. [TJ where 
the system is divided into three regions, A, B and C. 
We call all edges that go from one region to the other 
the boundary-crossing edges and label the subset of all 
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FIG. 1: (Color online) Example of a mathematical graph as- 
sociated to a physical graph state. We have displayed a possi- 
ble partition of this graph, splitting the system in three parts 
A, £>, and C. The vertices and edges in grey corresponds to 
the boundary qubits and the boundary- crossing edges respec- 
tively. 



such edges by X . All qubits connected by the boundary- 
crossing edges are in turn called the boundary qubits and 
the subset composed of all of these is called y. 

Open-system dynamics. - Our ultimate goal is to quan- 
tify the entanglement in any partition of arbitrary graph 
states undergoing a generic physical process during a 
time interval t. The action of such process on an ini- 
tial density operator p can be described by a completely- 
positive trace-preserving map A as pt = A(p), where 
p t is the evolved density matrix after time t. All such 
maps can be expressed in a Kraus representation, A(p) — 
"Yli^P^K^pK^, where y/pJ^K^ are called the Kraus opera- 
tors (each of which appearing with probability which 
satisfy the normalization conditions TrfiTTj^] = \ &nc \ 
Y^pPfi = 1 |13j . The Kraus representation guarantees 
that the map is (completely) positive and preserves trace 
normalization. When the map can be factorized as the 
composition of individual maps acting independently on 
each qubit, the noise is said to be individual (or indepen- 
dent); if not, it is said to be collective. 

A very important class of processes is described by 
the Pauli maps, separable (non-entangling) maps whose 
Kraus operators are given by tensor products of Pauli 
operators X, Y, Z, and the identity. Examples of these 
are the collective or individual depolarizing, dephasing 
or bit-flip channels |13| . As we show next, it is possible 
to determine the exact entanglement evolution of graph 
and graph-diagonal states (whose formal definition is pro- 
vided below) subject to individual Pauli maps. 

Exact entanglement of graph states under Pauli maps.- 

Let us start by recalling that a graph state is the si- 
multaneous eigenvector - of eigenvalue 1 - of the N gen- 
erators of the stabilizer group, that is, of the N operators 
consisting each of which of one X acting on each single 
qubit and Z's on all its neighboring ones pQ. Therefore, 
the application of an X or Y operator on a qubit A; of a 
graph state is equivalent to the application of Z opera- 
tors on all neighboring qubits of k, or on all of its neigh- 
boring qubits and on k itself, respectively. The action 
of any Pauli map A on a graph state is thus equivalent 
to that of another separable map, A, whose Kraus op- 



erators are obtained from replacing in the latter 
each X and Y operators by tensor products of Z and 
identity operators according to the rule just described 
|15j . Thus we need to consider how a general combina- 
tion of Z operators acts on a graph state. We use the 
multi-index p — (pi, pn), with pi — {0, 1}, to denote 
such a combination through Z^ 1 ® Z^ <g ... <g) Z^ N . The 
action of such operator on a graph state |G(y,£) ) gen- 
erates another graph state |G(y,£)-), orthogonal to the 
former one [TJ 110] . These considerations imply that pt 
can be expressed as 

Pt = A(|G (V)f)o )) = A(|G (v , 5)o » 

= £p*l G <v^K G cv,*) A |. (2) 

All possible 2 N graph states |G(y,£)-) associated to the 
graph G(y £) form a complete orthonormal basis of the 
iV-qubit Hilbert space. State ([2]| is a graph-diagonal 
state. Calculating the exact entanglement in any par- 
tition of the such state is in general a problem that in- 
volves an optimization over the entire parameter space 
of p t . In what follows we will show that it is possible to 
greatly reduce the complexity of this optimization prob- 
lem. Consider any partition of the state pt- We now 
factor out explicitly all the CZ gates but those corre- 
sponding to the boundary-crossing edges and write the 
state as 

Pt = CZ ij y ^2 l P 1 A G {y,x) 1 ){G(y,x) 1 \ 

{iJ}<££/X 7 ,<5 

® \g(v/y) s ){9(v/y) 5 \ CZ H , (3) 

{k,l}e£/X 

Here we have grouped together all indices inside p into 
two new multiple indices, 7 and 5. Multiple index 7 ac- 
counts for all possible graph states \G(y^x) 7 ) generated 
by applying tensor products of Z and identity operators 
to the graph state \Gry, X ) Q ) = ® { . ho}eX CZ i} ® \gry) Q ), 
associated to the boundary graph G(y^x) — {y, with 
\g(y) ) = l+i)- Multiple index 8 on the other hand 

accounts for all states \g(v/y)g} generated from Z or iden- 
tity operators on the state |<7(v/^) ) = ®iev/y °^ 
the non-boundary qubits V/y. Probability p 7j( 5 is defined 
as the sum of all p M such that K^\G(y iX ) ) ® \9(v/y) ) — 
\G(y,x) ) ® \9(v/y) s )- Because the CZ gates explicitly 
factored out in state ^ are local unitary operations with 
respect to the partition of interest, the entanglement of 
p t , E(p t ), reads 

EiJ^Pi^Giy^^iGiy^)^ ® \g(v/y) s )(9(v/y) s \)- (4) 

7,5 

where E is any convex entanglement quantifier not in- 
creasing under LOCC. In what follows, we first establish 
a lower and upper bound to this expression and then 
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show that these bounds coincide, obtaining the exact ex- 
pression of the graph-state entanglement evolution. 

First, consider an LOCC protocol consisting of measur- 
ing all the non-boundary qubits V/y of the state within 
brackets in Eq. Q in the product basis composed by all 
orthonormal states {\grv/y) s )} an d tracing out the mea- 
sured subsystem after communicating the outcomes. The 
remaining subsystem y is flagged by each measurement 
outcome S - meaning that outcome S provides full infor- 
mation about to which state y has been projected after 
each measurement run. The final entanglement after the 
entire protocol is then given by the average entanglement 
over all measurement runs. Since E is non-increasing un- 
der LOCC, E(p t ) must satisfy 

E(p t ) > J2psE(^2Ph\s)\G(y, Xh )(G {yXh \), (5) 

S 7 

where ps = ^2 7 P-y,8 is the total probability of occurrence 
of an event S and Pms) is the conditional probability of 
an event 7 given that event S has happened. 

On the other hand, convexity of E implies that 
E(p t ), as given by Q, must necessarily be smaller 

or equal to J28Ps E (j2-yPh\5)\G(y,x) 7 )(G { y iX)i \ (8 

\9(v/y)$)(g(v/y) s \j > which, since locally added ancillary 
systems do not change the entanglement, is in turn equal 
to the right-hand side of This means that the right- 
hand side of Eq. ^ provides at the same time an upper 
and a lower bound to E(p t ) and therefore yields its exact 
value, i.e. : 

E{ P t) = Y.Pt E (Y.PW)\G(y, Xh ){G ( y, X)l \). (6) 

s 7 

A comment on the implications of this exact result 
on the computational-cost is now in place. The cal- 
culation of the entanglement of systems composed by 
N = Ny + N v iy qubits (being Ny and iVy /y the num- 
ber of boundary and non-boundary qubits respectively) 
is a problem that, in general, involves an optimization 
over 0{2 2N ) real parameters. Through Eq. ^ such cal- 
culation is reduced to that of the average entanglement 
over a sample of 2 Nv /v states (one for each measurement 
outcome S) of Ny qubits, which involves at most 2 Nv /v 
optimizations over 0{2 2Ny ) real parameters. Thus the 
present method provides an exponential decrease in the 
computational power needed to calculate E(p t ), since 
only the boundary qubits appear in the computation of 
Eq. (§. 

In order to illustrate the power of the method we have 
calculated, using Q6J) , the exact entanglement of forma- 
tion Ep [16] of 1-D graph states under the action of in- 
dependent depolarizing channels, which mix, with prob- 
ability p, any one-qubit state with the maximally mixed 
state 1 /2 jT3 . In Fig. [2] we display the curves cor- 
responding to the bipartition first qubit versus the rest, 
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FIG. 2: (Color online.) Entanglement of formation (Ef) in 
the partition of the first particle versus the rest for 1-D graph 
states of 2 (black), 4 (grey) and 7 (red) particles undergoing 
individual depolarization as a function of the depolarization 
probability p. The dashed curve represents a size-independent 
lower bound. 



although other partitions can be considered. The 1-D 
graph state (also called the linear cluster state), given 
by \LC) = (g)^ 1 CZ iti+1 (g)f \+ k ), evolves from p = 
towards a final maximally mixed state at p = 1. Not 
only this calculation would have been impossible had we 
attempted a brute- force optimization approach, but also, 
since in this particular case the boundary qubits are just 
two, the use of (JsJ) allows to perform the calculation with 
no optimization at all, for an explicit formula for the en- 
tanglement of formation exists for arbitrary two-qubit 
systems [IB]- 

Beyond graph states and Pauli maps. - The expression 
^ is actually a method for calculating the entanglement 
of any graph-diagonal state as the one in ^ . Since Pauli 
maps acting on initial graph-diagonal states also produce 
graph-diagonal states, all the arguments used so far are 
also valid for this class of initial states. Furthermore, any 
quantum state can be depolarized to a graph-diagonal 
state by means of LOCC [TJ. Using again the fact that 
the entanglement of a state does not increase if an LOCC 
protocol is applied, one can see that the present method 
also provides (in general non-tight) lower bounds to the 
decay of the entanglement of any initial state subject to 
any decoherence process. 

Robustness of graph-state entanglement.- The devel- 
oped techniques can be further simplified to obtain new 
lower bounds to graph-state entanglement during all the 
evolution that, despite not being tight, can be calculated 
in a much more efficient way than §5§ and often turn out 
to be independent of the total number of qubits. This 
dramatically simplifies the study of the entanglement ro- 
bustness of graph states as a function of the system's size, 
a central question for the applicability of these states as 
quantum information resources. As an illustration, we 
compare next graph states of different sizes under the ac- 
tion of general 7V-qubit Pauli maps A that scale with N 
in a way such that, for each ^Jp^K^, the Kraus operators 
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of the map acting on M more qubits are obtained as ten- 
sor products of yJp^K^ with Pauli or identity operators 
on the other M qubits, weighted with some new proba- 
bilities that sum up to one (for each fi). That is, so that 
the total probability of event fi on the N first qubits, p^, 
remains the same. All the individual or collective Pauli 
maps mentioned above fall into this category. The state 
between brakets in Eq. Q can then be written also as 
Y,- ( P~f\ G (y,x) 7 )(G(y,x) 7 \ ® T l sP(sh)\9(v/y) s ){9(v/y) 5 \), 
where p 1 = J^gPi, s an d PtsM is the conditional proba- 
bility of 5 given 7. By tracing out the state of the non- 
boundary qubits (i.e., by disregarding the flag that lead 
to ([5]) above) and using again the fact that E does not 
increase under LOCC, we arrive at 

E{ Pt ) > E(J2P^\G(y, Xh )(G ( y !Xh \). (7) 
7 

Now, notice that - for the maps here-considered - proba- 
bility p 7 depends only on the boundary graph Gry^x) and 
the number of non-boundary qubits directly connected 
to it (the boundary graph is affected by the noise on 
up to its first neighbors), not on the total system size 
N. Bound ([7]) is unaffected by the addition of M extra 
particles if these new particles are not connected to the 
boundary subsystem. In the latter sense, and for the con- 
sidered noise scenario, noisy graph-state entanglement is 
thus robust with respect to the variation of the system 
size provided Gry t x) an d its connectivity to the rest do 
not vary. 

Size-independent bound Q (for the case of E — Ep) is 
compared with the exact entanglement, again for a linear 
cluster and the individual depolarizing channel, in Fig. [2] 

Discussion.- To summarize, in this work we have pre- 
sented a general framework to study the entanglement 
decay of graph states under decoherence. It is impor- 
tant to emphasize that any function that satisfies the re- 
quirements of convexity and monotonicity under LOCC 
falls into the range of applicability of the machinery here- 
developed. This includes genuine multipartite entangle- 
ment quantifiers, as well as those functions aiming at 
quantifying the usefulness of quantum states for given 
quantum informational tasks. 

To conclude with, let us make the following observa- 
tions. First, the techniques developed to obtain perfect 
bounds can also be applied to tackle some cases other 
than Pauli maps. For example, for graph states in the 
presence of individual thermal baths at arbitrary tem- 
perature, an LOCC procedure similar to the one used to 
obtain the bound jHJ, but using general measurements 
instead of orthogonal ones, can be used to obtain highly 
non-trivial entanglement lower bounds. 

Second, bound ([7]), when restricted to bipartite entan- 
glement, provides the same type of lower bound as the 
one used in section V-B of Ref. [TO] to find lower bounds 
to the entanglement lifetime for the case of E being the 



negativity. The present bound has the advantage of deal- 
ing with other possible partitions and general entangle- 
ment quantifiers. All these topics will be touched upon 
elsewhere. 
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